Invariant variational principle 
for Hamiltonian mechanics 



Alexey V. Golovnev^, Alexander S. Ushakov^ 



1. Ludwig-Maximilians University, astroparticle physics and cosmology de- 
partment. Theresienstr. 37, D-80333, Munich, Germany (on leave from 
St. Petersburg State University). Alexey.Golovnev@physik.uni-muenchen.de 

2. St. Petersburg State University, V.A.Fock Institute of 
Physics. Ulianovskaya ul. 1, Petrodvoretz, Saint-Petersburg, Russia. 
asushakov@gmail.com 



Abstract 

It is shown that the action for Hamiltonian equations of motion can be brought 
into invariant symplectic form. In other words, it can be formulated directly 
in terms of the symplectic structure uj without any need to choose some 1- 
form 7, such that uj = d'-f, which is not unique and does not even generally 
exist in a global sense. 



1 Introduction 



Hamiltonian mechanics is defined by a smooth function (Hamiltonian) on 
an even-dimensional manifold M^" equipped with a symplectic structure, i.e. 
a closed non-degenerate differential form. 



UJ 



uj^u{x)dx^ A fix'' = i ^ u^udx^ A dx" , x G M^''. 



tJ.>U 



li,u 



Using Darboux transformation these coordinates may be devided into two 
groups, X = (g^, q"',pi, ■..,Pn) with g* and Pi being generalized coordinates 
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and momenta respectively, but this transformation is not unique and in gen- 
eral can be performed only locally. 

Equations of motion are given by Poisson brackets 

2n 

= {x, H} = Y^ Lj^"'{x)d^H, (1) 
i/=i 

where the right-hand side can be regarded as a vector obtained by contrac- 
tion of Poisson bivector (dual to u) with the differential form dH; in local 
coordinates the bivector is represented by u^", the matrix inverse to Ufj^^. In 
the case of the standard symplectic structure, i.e. 

jiu I On En \ 

~\-En On J ' 

where En and On are the unit and the zero matrices respectively, a simple 
action exists for these equations of motion [T]: 




with the boundary conditions q{ti) = qi, 5(^2) = ^2 ■ 

But the expression does not seem satisfactory from geometric point 
of view because it is suitable only in Darboux coordinates and does not use 
the symplectic form explicitly. So, our goal would be to find an invariant 
geometric form of the action. In Sections 2, 3 we discuss dynamical systems 
with exact symplectic forms and, first of all, the very important case 
of "classical" Hamiltonian systems (i.e. those in M^" with the standard 
symplectic structure). The special emphasis is placed on the problem of 
boundary conditions [2]. After that we explore non-exact symplectic forms 
(see also [3] for Kahler manifolds as phase spaces, which are unitary phase 
spaces in the terminology of [3]). To the best of our knowledge there is no 
invariant variational principle in the literature for the general case. However, 
such systems may be of interest for different reasons. For example, systems 
with gyroscopic forces can not be presented in a straight-forward way as 
Hamiltonian systems with one- valued Hamiltonians The problem can 
be solved [4, 5] by some non-exact symplectic structure. In particular, the 
so-called Kirchhoff type systems [5] can be brought into Hamiltonian form 
on a cotangent bundle over 2-sphere with a volume form of the configuration 
space iS^ being added to the standard symplectic structure. It is enough to 
mention that systems of this kind include [5] a rigid body in an ideal incom- 
pressible fluid at rest at infinity, a rigid body with a fixed point in axially 
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symmetric potential field, spin dynamics in the A-phase of the superfiuid '^He 
(Leggett equations), etc; and at the quantum level we would have a sphere 
with non- commutative coordinates. Motivated by these reasons, in Sections 
4, 5 we establish an invariant variational principle for arbitrary Hamilto- 
nian systems. The last Section 6 is devoted to some examples and discussions. 



2 Brief review of known results 



As it was mentioned above, for classical Hamiltonian systems one can 
use the non-invariant variational principle with the action ([2]) to obtain the 
equations of motion ([T]). This action can be brought into invariant form as 
follows [HlEj: 

S = (^^u^, x^x^ -H^dt = j^' {p,q' - p,q') - dt, (3) 

= 1, . . . , 2?T,, i = 1, . . . ,n, where the summation over repeated indices 
is assumed and Dirichlet boundary conditions {6q{ti) = 5q{t2) = Sp(ti) = 
5p{t2) =0) on both coordinates and momenta variations are imposed, which 
is two times as much as the number of possible independent boundary con- 
ditions on the values of coordinates and momenta themselves in Cauchy or 
boundary problem for the equations of motion. We have a mismatch between 
the number of required variational boundary conditions and the number of 
independent solutions of equations of motion. 

So, the Hamiltonian variational principle differs in this sense seriously 
from the Lagrangian one. Actually, one may suspect that it can be formulated 
in a more usual way by fixing the values of the functions cbi = arctan — at 

Pi 

the boundary points. Indeed, variations of 0-s equal 

PiSqi-qi6pi 

d(pi = — — 2 — ; 

QI + PI 

and S(j)i = at t = to and t = ti is exactly what we need to make the 
boundary terms in the first variation of equal zero. Equations (pi = const 
define some n-dimensional submanifolds in 2n-dimensional phase space. And 
it would be curious to introduce new pairs of coordinates in the following 
way: 

qi = Ti sin (pi Pi = Ti cos (pi. 
With the definition Pi = rf /2 it is easy to get ^{piqi — qiPi) = Pi(pi and 

-12 



S = ^J2P^<i'^-H{<P,P)]dt. 



It looks quite good. But the point is that these variables are nothing more 
than a particular choice of Darboux coordinates. We could well have started 
with the action (IHl) for these canonical variables and considered other prefered 
functions of the form 0j = arctan ^. Hence we should stress once more that 
there is no invariant way to fix only a half of boundary conditions without 
sticking to some coordinate choice, because what we have done here is just 
a canonical transformation to new variables, Pi and 0^. 

It is of crucial importance for the geometric meaning of the action ([3]) 
that the phase space is M^" because the radius-vectors are used. (But the 
symplectic form could be arbitrary.) However, there is a natural analogue of 
the action ([2]) which is valid for any exact symplectic manifold. (Note still 
that for any compact symplectic manifold (without boundary) uj is always 
not exact [H].) 

Suppose that the symplectic form is exact: u = d^y, where 7 is some 1- 
form (for classical systems 7 = Pidq^). (Even if the symplectic form is not 
globally exact, such 1-forms always exist locally.) Then the matrix u^i, can 
be expressed in terms of the components of 7 = A^dx^: u^y = d^A^ — d^A^. 
And the action can be written as follows [2]: 



Actually, the action (jll) is just a special case of Birkhoff theory (Pfaff-Birkhoff 
variational principle) [3], [9], HD] 



with B staying for Hamiltonian. The equations of motion should be derived 
from the variational principle with Dirichlet boundary conditions. This for- 
malism is valid also for non-conservative systems with time-dependent sym- 
plectic structure. 

The variation of the action (jlj) yields Hamiltonain equations of motion if 
one imposes Dirichlet boundary conditions as it was done above in R^" with 
the same problems persisting. In principle, provided that we know an explicit 
form of A^ we can introduce an analogue of the functions 0. But if we want 
to allow for different possible choices of 7 then we are forced to admit the 
whole set of boundary conditions, 5q{ti)=5q{t2)=5p{ti)=5p{t2)=Q. 

This mismatch between the number of boundary conditions for the vari- 
ational problem and the number of boundary conditions for the equations of 
motion often makes physicists to worry about more that it really deserves. 
The habit to identify one kind of boundary conditions with another is so 
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strong that many distinguished authors fHX [T^ [TB] do not fix initial and 
final points for the action ([2]) completely but impose only one half of the 
conditions {Sq(ti) = 5q{t2) = 0), even if they are going to discuss canonical 
transformations in the same text. 

Actually, the story of this approach goes back to the early days of 
Hamiltonian mechanics. But already in the book by Poincare [13] a dis- 
cussion of canonical transformations appeared with the conclusion that 
such principle has different forms in different canonical coordinates (in the 
sense of different boundary conditions). For this reason it is commonly 
accepted in mathematical literature to fix the boundary points completely 
by 4n conditions [15j, thus providing the variational principle with a kind 
of invariance. And it is pleasant to note that the same is done in the 
classical textbook for physicists [TB] (and in some other physical books 
[T7] ) where it is also stressed that only the whole set of boudary conditions 
allows one to add a total time derivative term ^/(p, q) to the integrand in ([2]). 

3 The case of "wrong" boundary conditions 

In this Section we would like to add several comments and to give a 
variational characteristic of correctly chosen boundary points, which is in a 
sense obvious and presumably not easy to use but, may be, worth mentioning 
still. Some other problems would be discussed in a forthcoming article by L.V. 
Prokhorov and A.S. Ushakov. 

In any case the correct statement is that the physical trajectory brings the 
action to the stationary value in the class of trajectories with fixed bound- 
ary points in the phase space. But an important fact about the variational 
principles is that one would probably wish to use the relevant principle in 
order to get the equations of motion or to apply it instead of the equations if 
it can make the mathematics easier. At this point we encounter with a kind 
of problem: the number of boundary conditions required is greater than the 
number of Hamiltonian equations of motion. And for every initial point at ti 
there is only one final point at ^2 for which the extremizing trajectory exists. 
To get the physical trajectory one needs to guess the unique final point from 
the continuum of all the phase space points. But this problem can't lead to 
a false trajectory because if the final point is chosen wrong the equations of 
motion can not be satisfied and the action has no critical point in our class 
of variations. One could even say that it is not a problem at all because one 
can choose arbitrary initial point and fix the final point formally without any 
idea of where it is fixed [2]; it will allow him to get the equations of motion 
and find the location of the final point after that. 
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The well-known and very important property of the action ([2]) is that 
it is bounded neither from below nor from above. As a consequence, the 
physical trajectory is not a kind of extremum but rather is a saddle point. 
In particular, it means that the direct methods of calculus of variations can 
not be directly applied. And it was not until about 30 years ago that certain 
progress has been made in this direction [THl [IHl El 1201 due to technique 
of approximation by some finite-dimensional critical point problem [TH]. It 
means also that if one shifts a final point a little bit from its correct position 
then the range of the action does not change radically, just its saddle point 
is gone. 

We can regard the first variation of the action ([2]) as a linear opera- 
tor acting from L^-space of functions 6x{t) : [^1,^2] ~^ I^^" to M, which 
depends on a chosen path from x{ti) to x{t2). Its operator norm equals 

\\5s\\ = JfdtEi ((^^ - Iff)' + + If )')• o^ly physical 



trajectory we have \ \Ss\ \ = 0, otherwise H^^H > 0. If we have a "wrong" final 
point then for all the possible paths \ \^s\ \ > 0. Moreover, the minimal value 
of ll^sll normally exists, so that it is separated from 0. 
The minimal path for the functional 1 1 ^5 1 1 is given by 

d f dH\ d^H ( dH\ d^H ( dH\ 

d f dH\ d'^H ( dH\ d^H f dH 



dt \ dqi J dqidpk \ dqu J dpidpk \ dpk, 

which allows exactly An independent boundary conditions. If these conditions 
are "right" we have also g.j — |^ = 0, + |^ = and H^^H = 0. In other 

cases 1 1^5 1 1 > 0. Let us consider a very simple example with n = 1 and 
H = implies H^s'lP = / dt {{q — pY + p^) . The minimum is given by 



q{t) = + + (C3 - Ci)t + C4 and p{t) = + Cat + C3. For the 
initial point {qo,po 7^ 0) we have to choose (go +^0(^2 — ti),Po) as the final 
point. But if we take (go + 0^0(^2 — ^1)5^0 + P): the constants Ci would be 
C, = go, C, = po, C, = ,,l^t?:mXv ^2 = lAr - finally 
\\5s\? = Cl {^{h - hf + ih - tl)) + Only for « - 1 = /? = do we 

have \ \6s\ \ = 0. 

Of course, we used here even more complicated equations of motion than 
the original ones. But the variational problem is quite different from that of 
([2]) because | l^s] | has a unique global minimum, and its value in principle can 
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Fig.l ^ 



be found by direct methods of variational calculus. After that, these minimal 
values can be considered as a function of the final point (the initial point is 
fixed) which, in turn, also has a unique global minimum in the "right" final 
point. The minimal value is zero. 

Note that if the symplectic manifold is different from M^" we may 
consider the variations as elements of the tangent space. We don't want to 
go into any details here but for any symplectic manifold there exist (not 
unique) an almost complex structure and a Riemannian metric which are 
in a sense compatible with oj (see, for example, [20], p. 14), and it makes a 
room for the constructions explained above to be used in the general setting 
of the action (jll) . 



4 Invariant variational principle 

In the general case one should use the symplectic form uj explicitly. And 
all the vectors we have are in the tangent space. It means that the surface of 
integration should be 2-dimensional. Looking at the formulae ([3]) and (jlj) we 
expect that the action should have the following form: 

S = j {u-dHAdt) = J (^^tOf.^dx" A dx" -dHA dt^ . (5) 

More precisely, we consider a trivial bundle over M^" : _p2n+i _ ]\^2n ^ ]^ 
(extended phase space) and endow it with a new coordinate t (time) so that 
the basis of 1-forms gains one more, (2n+l)-th, element, dt, with an axiom 
dt Adx = —dx A dt. The differential forms u and dH are defined in F simply 
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by their coordinate expressions in M which may be invariantly interpreted 
as a pullback of these forms generated by the natural projection of F to 
M. Any one-parameter family of initial points e G [0, -E] defines a 

two-dimensional surface a of trajectories in F with the curve of the final 
points at t = ^2- It is important that the initial and final curves (71 and 
72) should be transversal to the physical trajectories for the surface a to be 
well-defined. One possible way to ensure it is to choose the curves transversal 
to the hypersurfaces of constant Hamiltonian. In this case the parameter e 
gains also a possible interpretation as a value of the Hamiltonian function, so 
that one considers a family of trajectories with different amounts of energy. 
In the rest of the Section we prove under these assumptions the following 
Theorem: 

1. The surface a is a stationary one for the action ^ in the class of 
smooth surf aces X = x{e,t), e G [Oj-E], t G [^1,^2] with fixed ends 6x{e,ti) = 
Sx{€,t2) = 0, Ve (Fig. 1). 

2. An arbitrary surface a from this class is stationary if and only if its 
boundary trajectories x{0,t) and x{E,t) are physical, i.e. they satisfy the 
equations of motion. 

Let us start with a simple case when the surface a and small variations 
of it belong altogether to one coordinate chart of the phase space manifold. 
Strictly speaking, it means that we use only a contractible domain in M^" 
and there does exist some suitable 1-form 7. From its definition = (i7 it 
follows that the action ^ equals to the action (jH) on da and the statement 
is trivial. Nevertheless we want to proceed with explicit calculations in order 
to make an illustration of how it works if one does not know an appropriate 
1-form 7. We perform the variations here in a somewhat formal way. We 
vary coordinates even under the differential symbols in differential forms 
as if they were just ordinary functions. It can be justified if we take into 
account that for the small variations of the surface there is a natural one- 
to-one correspondence between points of the initial and the final surfaces. It 
allows us to vary the integrand instead of the domain of integration. It is 
also important to mention that we consider only continuous variations of the 
surface a so that all the surfaces are homotopic to each other and the Stokes' 
theorem is applicable. 

The variation of the action (I5l) is equal to 

5S= f {^d^ujf.M'^dx^ ^dx'' + ]-uj^^da{5x^)dx'' ^dx^+ 

+ ]^uj^udc{6x'') dx^" A c/x° - d,3{daHSx")dx'^ A dt^ . 
We notice that the second and third terms under the integral are equal, and 
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integration by parts gives 

+ {daiuj^iydx'') - daUJ^„6x'^)dx" A dx"- 

- d(i{daH6x'')dx^ A dt 

Using the Jacobi identify 

we get 

6S = I (^df.UauSx^'dx^ A dx" - ^d^uJafiSx^'dx^ A dx''+ 

Jcr ^ 

+ ipai^y^ybx^) - daUJ^,y6x'')dx" A dx"- 
- dfs{daH6x'')dx^ A d?j = 

= I ( ~ -d^uJauSx'^dx^ A dx" + -duUJa^6x°'dx'' A dx^+ 

+ d^{u^Jx^')dx'' A dx"" - dp{d^HSx'')dx'^ A dt 



^dc,{u^Jx'')dx'' A rfx^ - d^{dpH6x^)dx'^ A rft 



uj^jydx^dx" - dpEdx^^dt 



9{[0,£;]x[ti,t2]) 

t2 



uo^ydx^dx" - dfj_H6x^dt 



j ' - d^H{x{E,t))\5x^dt- 



t2 
tl 

The variations 5a;'^ are arbitrary and can be performed independently at 
e = and at e = ii^. Hence, 55* = is equivalent to Hamiltonian equations of 
motion uj^j^x'^ = d^H for x(0, t) and x{E, t). It means that in the space of all 
possible one-parameter families of curves we have an "infinite-dimensional 
stationary manifold" for the action (I5l). The intermediate trajectories do not 
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influence the stationarity condition because the differential form under the 
integral in is closed and any variation within one homotopy class with 
the whole boundary da fixed is subject to the Stokes' theorem. 

It's worth to mention that before the last step of the calculation all x and 
6x could be regarded as functions of point in F^n+i^ Q^jy ^^j^g jg^g^ g^gp 
we restrict ourselves to the boundary lines and consider these functions as 
well-defined functions of time. It means that we could vary the time in 
independently. It is easy to see that the consequence of such variation is that 
the Hamiltonian does not change with time (along the physical trajectory). 
It follows also from the equations of motion. It's not surprising that we have 
got nothing new because the variation of time is equivalent to some variation 
of dynamical functions x{t). 

Suppose now that the surface a in our calculation is divided into several 
parts by a number of internal lines which are transversal to each other and to 
da. For each part the variation of the action yields the equations of motion on 
its boundary. Now we want to glue these parts together. How can the result 
for initial a be restored after that? Clearly we have to add up the integrals 
over all the parts and demand that the variations of different parts should 
be equal each other on those boundaries which are going to be identified. 
Then integrals over all internal lines in 6S come up twice with opposite 
signs and cancel each other. (We know that it really had to be so due to 
the Stokes' theorem.) Again we get the Hamiltonian equations of motion for 
the boundary trajectories of a. It allows us to prove the Theorem in its full 
generality. Indeed, even if we can't embed a surface a into a contractible 
domain in M^", we still can divide it into small parts such that every part 
together with its nearest neighbours belongs to some contractible domain. 
The action functional is invariant under coordinate transformations. Each 
part of the surface can be varied in any coordinate system and contributions 
of division lines cancel each other. The final result contains only equations 
of motion on the boundary trajectories of a. The theorem is proven. 

After this analytic derivation we would like to add a nice geometric picture 
of the theorem. As it was already mentioned, those variations which do not 
change boundary trajectories are irrelevant due to the Stokes' theorem. Let 
us consider a variation of the surface a which is not trivial only in a small 
vicinity of a part of one boundary trajectory. Then the variation is just the 
difference between two integrals over small pieces of a and a = a + 6a. And 
due to the Stokes' theorem again this difference is equal to the integral over 
some surface which connects the initial trajectory and its image after the 
variation. So, it means that the integral of u — dH A dt over arbitrary thin 
stripe along the physical boundary trajectory should be of order 0{{5xY) 
with no regard to the orientation of the stripe. (Note also that we can take a 
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Fig. 2 



not too long part of the trajectory and use the Darboux coordinates on the 
stripe if we hke.) It means that [uj — dH A dt){l, a) = for any vector a and 
a vector / tangent to the trajectory. It is precisely the equations of motion. 

Let us also briefly mention that we could define a "superextended " 
phase space with two new coordinates, H and t. In this case we should 
perform variations only on a hypersurface H = H{x) and may consider the 
parameter e as a coordinate along the H-axis (Fig. 2). 



5 Some additional remarks 



Curiously enough, the idea of some 2-form integration over 2-dimensional 
surface between two trajectories appeared recently in [21] which goes in a 
direction somewhat opposite to the lines of our work. Systems with velocity 
dependent forces, generally admiting neither Hamiltonian nor Lagrangian 
formulation (see, however, [22] and [23]), are considered in this reference 
in the language of positions and velocities (we would like to remind here 
that, on the contrary, our main goal is the treatment of Hamiltonian systems 
for which it is generally impossible to separate coordinates and momenta). 
For these (generally dissipative) systems a variational principle is obtained 
which yields the equations of motion and some more equation on the bulk of 
the integration surface with unclear dynamical meaning. So that a difficult 
problem of joint solvability arises, see [21] for details. In spite of certain 
similarity, the intersection of this principle with ours is rather trivial. In 
the case of conservative systems it reduces to the Lagrangian version of the 
action ([5]), but the relevant systems are those for which the coordinates and 
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momenta are clearly separated, and much simpler action principle of the form 
(12]) or (HD exists. 

Note also that we could formulate the variational principle invariantly 
but without dealing too much with the language of exterior calculus. The 
obvious relation dx'^ = x^dt + x'^de leads us to the action integral: 



And varying it with respect to x^{t,t) in the same way as in the Section 4, 
we get the same result 



Then we have to mention the problem of boundary conditions again. It 
may seem to be even more intricate in the case of non-exact forms because 
if the initial points in the phase space are given one has to guess properly 
the final points at least for two trajectories (at e = and e = E) for the 
stationary surfaces to exist. But we can easily reformulate our principle: 
choose in _F2ra+i Qj^jy physical trajectory and one auxiliary line with the 
same initial and final points such that a nonsingular surface in the phase 
space exists with the boundary equal to these two curves. Then we can take 
our action integral along the surfaces of that kind with the auxiliary 
line fixed and the physical trajectory free to change (except the boundary 
points, of course). By literally the same calculations as in the Section 4 it 
can be easily verified that all surfaces for which the equations of motion are 
valid on the trajectory would compose the stationary manifold for the action 
considered. 

And a final remark deals with Maupertius principle. Suppose we consider 
only the first term in the action ([5]) but vary this action only in the class 
of surfaces with boundary trajectories laying on two different hypersurfaces 
of constant Hamiltonian. The result is that for every vector l'^ tangent 
to the H = const hypersurface in M^" the equation uj^^^x^l" = should 
be held true on the boundary trajectories. In the phase space M^" this 
equation uniquely defines a curve which is compatible with the equations 
of motion a;(x, .) = —dH{.) because dH{l) = for the vectors considered. 
But the time coordinate may be chosen arbitrarily because the equation is 
reparametrization invariant. So we get the invariant form of the Maupertius 
principle. 





e=E 
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6 Examples and discussion 



The usual action principle for Hamiltonian systems is given by (jlj) (or by 
([2]) for classical systems). On exact manifolds the 1-form 7 is globally well- 
defined and the action (jlj) can be perfectly used, but one should remember 
that although the choice of 7 is a coordinate-free procedure, it is not unique 
and effectively means a kind of distinction between what we would call mo- 
menta and what we would call coordinates. For non-exact manifolds 7 exists 
only locally, and even if with a particular choice of 7 we were able to pursue 
a variational approach for one trajectory it doesn't mean yet that we would 
be able to repeat it for some another path without changing the set-up. 

Our variational principle is completely invariant and can be used for any 
Hamiltonian system. Let us consider a few very simple examples. First of all, 
we take a sphere with oj = cos 6* d9 Adip in spherical coordinates, 9 G [— |, |], 
ifi G [0,27r) and Hamiltonian H = sin6. (One can check that u is non- 
degenerate and H is smooth near the poles 6 = ±|, for example, by going to 
new coordinates x = cos 6 ■ cos ip, y = cos 6 ■ sin (p in which u = 

y/l-x'^-y^ 

H = — — ?/2 are obviously well-defined near x = y = 0.) The equations 
of motion are 9 = 0, (p = 1. Let ip{e, ti) = 0, 9{e, ti) = e, < e < _E < | be a 
family of initial data and (p{e,t2) = a{t2 — ti) , 6'(e, ^2) = e - a family of final 
points. (We use a to show what happens if final points are "wrong".) For this 
case (one degree of freedom) irrelevance of intermediate trajectories is obvious 
because the surface of integration x{e,t) is completely defined by boundary 
paths ^(0(t)) at e = and e = ip{E,t), (p{0,t), 9{E,t) = f2i(p{E,t)), 
9{0,t) = fi{(p{0,t)). The action after one half of the integrations have 
been performed, turns into the following form: 

Q{t2-ti) 




dip (sin(/2(¥?)) - sin(/i(^))) - 





- j dt{sm{f2{^{E,t)))-sm{f^{y,{0,tm. 

ti 

The variations of f2{'p),ip{E,t) and fi{(p),(p{0,t) should be performed in- 
dependently, and actually we can even fix one of the boundary trajectories 
and obtain the equations of motion only for the second path. If we vary 
ip{E,t) as a function of time we have 6^,3 = — f dt cos{f2{(p)) f2 S(p. It im- 
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plies /2 = const, i.e. 9 = 0. Now we vary f2{'^) independently: 



a{t2-ti) t2 

5f^S= J d(p cos(/2(<^))5/2(<^) - y dt cos(/2(v2))5/2(<^) 

tl 

and after changing the variable in the first integral according to dip = Lpdt we 
get if = 1. The consistency condition (Ayj = j (pdt) leads to a = 1. There is 
no stationary point for the action S otherwise. 

Actually, is a very simple manifold and one could use 1-form 7 = 
sin 9 dip on it. But this 7 is singular in the poles, and if we were not so clever 
we could take 7 with poles on the physical trajectory. One more remark is 
that we could make E —>■ ^ and convert the action (jl]) to a rather nice form 

S = J CO - j Hdt (6) 

where t2 — ti = 2tt and V is the upper part of the sphere with a periodic 
trajectory as a boundary. This expression is valid only for contractible closed 
trajectories and is being successfully used for the purposes of symplectic 
topology [20l I2I] , but it is not well known among physicists. (The Authors 
were unaware of it before searching the literature for historical references for 
this article.) We have to note that the actions and present even more 
intricate problems for the direct variational methods because if one changes 
the homotopy class of a or D he will not affect the stationarity of the action 
with respect to small continuous variations of the surface but will shift the 
value of the action by some element of the so-called period group which 
consists of all the real numbers obtained by integrating u over submanifolds 
homeomorphic to S^. Sometimes this ambiguity is only of academic interest 
because, for example, for systems on it means only that we could force 
the disk T> to wrap several times around the whole phase space. (Note that 
in this case the values of action ([6]) would aquire an additional constant, a 
multiple of the total area. It will change nothing for the variations of 

but no 1-form 7 would be correctly defined.) But unfortunately it is not 
always so good, and sometimes the period group may even be everywhere 
dense in M, see |20], p. 228. 

Needless to say, our principle works equally well for both contractible 
and non-contractible loops. Consider now a torus with u = d9 A dip, 
9 e [0, 2tt), (f G [0, 2tt), H = sm9 and equations of motion 9 = 0, ip = cos 9. 
(This example is oversimplified not only due to the local existence of j = 9 dip 
but also because for the torus one can use a trick [6] of considering the action 
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for the periodic functions in M^. But still we do not want to complicate 
it here.) Our periodic boundary conditions 9i = 9f = e E [0,E], ipi = 0, 
(fif = 27Tn, n G Z imply 

27rn T 


Suppose that we fix the e = path and vary only the e = E trajectory. Vari- 
ation of f2{<^) yields /2 = const and after that varying ip{t) we get ip = cos 6 
with the consistency condition again Aip = J ipdt. The last equation means 
that trajectories satisfying the periodicity relation (p{t + T) = (p{t) + 27m 
exist at cos 6' = There are several different solutions to the equation 
SS = for large enough T. (Different solutions have different boundary tra- 
jectories, and each one of these solutions also has an infinite degeneracy due 
to possibility of wrapping a any number of times around the torus without 
changing the boundaries.) In the invariant setting we do not have to worry 
about the location of these stripes a, but if we decided to use j = 9 dip it 
would be necessary to cut the torus along the line of 6^ = 0. 

These n = 1 examples do not allow us to illustrate the Maupertius prin- 
ciple. For this purpose let's consider the simplest example of n = 2 system: 
j^2n ^ ^2 ^ ]^2^ ^ ^ pQg^ de Adip + dpA dq, H = sinO + ^. The equations 

of motion are 6 = 0, p = 0, (p = 1, q = p. To apply the Maupertius principle 
we first find the surface of constant Hamiltonian: sin 6' + ^ = const. After 
that we take the boundary conditions: ipi = 0, ipj = t2 — ti, 9i = 9f = ae, 
qi = 0, qf = [3e{t2 — ti), Pi = Pj = /3e and the abbreviated action 

Sm= uj= / dip (sin(/2((/?)) - sin(/i(v?))) + dq (^2(9) - ^i(g)) 



where 9{E,t) = f2{ip{E,t)), p{E,t) = g2{q{E,t)) and similarly for fi,gi. 
Varying the second trajectory we have 6Sm = J dip cos(/2) 5/2 + J dq 5g2. 
We consider q = q{ip) and use the relation cos(/2) 5/2 + 02^92 = on the 
H = const surface. The result is 6Sm = J dip {—g2 + <i'{.v)) ^92 and q'{ip) = p. 
Independent variation of q{ip) yields / dip g25q' and g2 = const. Due to 
H = const it implies /2 = const and finally we have 6 = aE = const, 
p = j3E = const, ^ = PE = const. 

We have presented a new completely invariant approach to the variational 
formulation of Hamiltonian mechanics. Our principle can be applied to any 
Hamiltonian system, but for the direct methods of variational calculus it 
inherits all the usual mathematical difficulties related to action principles in 
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Hamiltonian form. Therefore the question of how far one could go further 
with it remains an open problem of mathematical nature. 

The Authors are grateful to anonymous referees for useful comments and 
suggestions. 

References 

[1] Arnold V.I., Mathematical methods of classical mechanics. Springer, 
New York, 1989. 

[2] Martinez-Merino A., Montesinos M., Annals Phys., 321, 318-330 (2006); 
preprint gr-qc/0601140 . 

[3] McEwan J., Found, of Phys., 23 (2), 313-327 (1993). 

[4] Souriau J.M., Structure des systemes dynamiques. Dunod, Paris, 1970. 

[5] Novikov S.P., Russ. Math. Surv., 1982, 37 (5), 1-56. 

[6] Conley C.C., Zehnder E., Invent. Math., 73, 33-49 (1983). 

[7] Sergi A., preprint Icond-mat/ObOSld^. 

[8] Bolsinov A.V., Fomenko A.T., Integrable Hamiltonian systems. 
Taylor & Francis Books Ltd, 2003. 

[9] Birkhoff G.D., Dynamical Systems. AMS, 1927, 1983. 

[10] Santilli M. Foundations of Theoretical Mechanics, I, II. Springer, Berlin, 
1978, 1983. 

[11] Landau L.D., Lifshitz E.M., Mechanics. Pergamon Press, 1976. 

[12] Konopinski E.J., Classical Descriptions of Motion. W.H. Freeman and 
Company, 1969. 

[13] Henneaux M., Teitelboim C, Quantization of Gauge Systems. Princeton, 
New Jersey, 1994. 

[14] Poincare H., Les Methodes Nouvelles de la Mecanique Celeste, V.3. Paris, 
Gauthier-Villars, 1899. Ch. 29. 

[15] Arnold V.I., Kozlov V.V., Neishtadt A.L, Mathematical Aspects of Clas- 
sical and Celestial Mechanics. Springer, 1997. 



16 



[16] Goldstein H., Poole Ch., Safko J., Classical Mechanics, 3-rd ed. Addison 
Wesley, 2002. Ch. 8, sect. 8.5 

[17] Fliessbach T. Mechanik. Spektrum Akademischer Verlag, 1999. K. 28. 

[18] Rabinowitz P.H., Comm. Pure Appl. Math., 31, 157-184 (1978). 

[19] Rabinowitz P.H., J. Diff. Eq., 33, 336-352 (1979). 

[20] Hofer H., Zehnder E., Symplectic Invariants and Hamiltonian Dynamics. 
Birkhauser, Basel, 1994. 



[21] Koclian D., preprint hep-th/0703073; see also Koclian D., 
|hep-th/0612115 for some geometric issues. 

[22] Gitman P.M. , Kupriyanov V.G., Eur. Phys. J. C, 50, 691-700 (2007); 
preprint hep-th /0605025[ 

[23] Gitman D.M., Kupriyanov V.G., J. Phys. A, 40, 10071-10081 (2007); 
preprint larXiv : 71 0.45321 

[24] Cieliebak K., Floer A., Hofer H., Math. Zett., 218, 103-122 (1995). 



17 



